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Abstract: Differential equations of infinite order are an increasingly important class of 
equations in theoretical physics. Such equations are ubiquitous in string field theory and 
have recently attracted considerable interest also from cosmologists. Though these equa- 
tions have been studied in the classical mathematical literature, it appears that the physics 
community is largely unaware of the relevant formalism. Of particular importance is the 
fate of the initial value problem. Under what circumstances do infinite order differen- 
tial equations possess a well-defined initial value problem and how many initial data are 
required? In this paper we study the initial value problem for infinite order differential 
equations in the mathematical framework of the formal operator calculus, with analytic 
initial data. This formalism allows us to handle simultaneously a wide array of different 
nonlocal equations within a single framework and also admits a transparent physical in- 
terpretation. We show that differential equations of infinite order do not generically admit 
infinitely many initial data. Rather, each pole of the propagator contributes two initial 
data to the final solution. Though it is possible to find differential equations of infinite 
order which admit well-defined initial value problem with only two initial data, neither 
the dynamical equations of p-adic string theory nor string field theory seem to belong to 
this class. However, both theories can be rendered ghost-free by suitable definition of the 
action of the formal pseudo-differential operator. This prescription restricts the theory to 
frequencies within some contour in the complex plane and hence may be thought of as a 
sort of ultra-violet cut-off. Our results place certain recent attempts to study inflation in 
the context of nonlocal field theories on a much firmer mathematical footing. 
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1. Introduction 



1.1 Differential Equations of Infinite Order in Theoretical Physics 

Differential equations containing an infinite number of derivatives (both time and space 
derivatives) are an increasingly important class of equations in theoretical physics. Nonlocal 
field theories with infinitely many powers of the d'Alembertian operator □ (given by □ = 
— 9| + in flat space) appear ubiquitous in string field theories (for a review 

see [1^]). This nonlocal structure is also shared by p-adic string theory |jl2| (see also 
fl^ ), a toy model of the bosonic string tachyon. Yet another example of a field theory 
containing infinitely many powers of □ can be obtained by quantizing strings on a random 
lattice 1 14] (see also |l^). Moreover, field theories containing infinitely many derivatives 

due to a wide 



have recently received considerable attention from cosmologists [16|- 
array of novel cosmological properties including the possibility of realizing quintessence 
with w < —1 within a sensible microscopic theory ||l^-[21|, improved ultra-violet (UV) 
behaviour [22, 23 1, bouncing solutions [^-|24] and self- inflation |25| (see also for 
applications both to bouncing cosmologies and also to dark energy and see ^ for a 
discussion of the construction of cosmological solutions in inflnite order theories). In [^] it 
was shown that cosmological models based on p-adic string theory can give rise to slow roll 
inflation even when the potential is extremely steep. This remarkable behaviour was found 
to be a rather general feature of nonlocal hill-top inflationary models in [^]. In [^4| it was 
shown that nonlocal hill-top inflation is among the very few classes of inflationary models 
which can give rise to a large nongaussian signature in the Cosmic Microwave Background 
(CMB). 

Differential equations with inflnitely many derivatives which arise frequently in the 
literature include the dynamical equation of p-adic string theory |12] 



P 



fl.l) 



where p is a prime number (though it appears that the theory can be sensibly continued 
to other values of p also |3^) and we have set rus = 1. A second popular example is the 
dynamical equation of the tachyon fleld in bosonic open string fleld theory (SFT) which 
can be cast in the form (see, for example, fsF 



[(1 + 



□ e 



where c = ln(3^/4^). In both cases the field (p is a tachyon representing the instability 
of some non BPS D-brane configuration. More generally, there is considerable interest in 



applications of a wide class of equations of the form [24, [2^, 27 

F (a) 4> = u {(t>) 



(1.3) 



where U{(j)) = V'{(j)) is the derivative of some potential energy function V{(j)) associated 
with the field 4> and we refer to F{z) as the kinetic function. Equations of the form ( pT^ ) 
are interesting in their own right from the mathematical perspective and some special cases 
have received attention recently [^ 
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Somewhat more general classes of infinite order differential equations, in which the 
derivatives do not necessarily appear in the combination □ = —d^ + V^, arise in noncom- 
mutative field theory fluid dynamics and quantum algebras [42|. 



1.2 The Importance of the Initial Value Problem 

Of particular interest is the fate of the initial value problem for infinite derivative theories. 



When does equation (1.3) admit a well-defined initial value problem - even formally, that 
is ignoring issues of convergence - and how many initial data are required to determine a 
solution? Such questions are fundamental to any physical application. To emphasize the 
nature of the problem that we are solving, let us consider a few trivial examples. First 
consider the equation 

{df - n?f = (1.4) 
It is easy to see that, due to the presence of the operator — m?, the function 

Mt) = Ae""* + Be-"^' (1.5) 

provides a two-parameter solution of ( |1.4| ). Of course, this is not the most general solution. 
To 4>i one could also append 

Mt) =t[Ce"'' + 06-""'] (1.6) 



Since equation (1.4) is fourth order in derivatives it is not surprising that the four parameter 
solution (f){t) = 01 (t) -|- 02(0 provides the most general possible solution and the initial 
value problem is well-posed with four initial conditions. We now consider the somewhat 
less trivial example 

Ai{-d?)M) = (1.7) 
where Ai(2;) is the Airy function. It is perhaps obvious that a solution is provided by 

M) = Yl K^""* + bne-"^"'] (1.8) 

n 

where {rn^} are the zeroes of the Ai{—z) (the first few roots are = 1.17, = 3.27, 
• • • ). What is probably less obvious is whether (^]^) provides the most general solution of 
(1.7). Could there be additional solutions which are not as trivial to construct as (|l.8|)? 



Since the original equation (1/7) was infinite order, it is not entirely obvious. An even less 



trivial example is provided by the equation 



52-m20 = O (1.9) 

which has appliciations to bosonization. How many initial conditions are necessary to 
specify a solution of (pTg])? 

In this paper we study the initial value problem for differential equations with infinitely 
many derivatives in the context of the formal operator calculus. The initial value problem 



for various nonlocal theories has also been studied in [44|-|5C]. In this paper we lay 



down a simple and intuitive formalism for studying the initial value problem in nonlocal 



-3- 



theories which is sufficiently general to handle a wide variety of different types of nonlocality. 
(In particular, our formalism is sufficiently general to handle the motivational examples 
( |1.4 |1.7| , 1.9).) Though we do not derive any new solution of pT^ ) we nevertheless feel 
that it is instructive to re-consider such equations in the context of our formalism. Our 
approach may be readily applied also to other nonlocal equations which arise in physical 
applications. 

It is not uncommon to see equations of the form (^]^) described as "a new class of 
equations in mathematical physics" in the string theory and cosmology literature. However, 
it happens that differential equations of infinite order have been studied in the mathematical 
literature for quite some time. Indeed, the study of linear differential equations of infinite 
order was the subject of an extensive treatise by H. T. Davis as early as 1936 |51]! This 



treatise and papers by Davis |52| and Carmichael [^3[ give an account of this theory as 
it stood at the time. The topic was further developed from an analytical perspective by 
Carleson [54|, who showed that the solutions of differential equations of infinite order need 
not be analytic functions, and obtained sufficient conditions in terms of the coefficients 
of the equation for the solutions of the initial value problem to be analytic. We should 
also mention that initial value problems for some special classes of differential equations of 
infinite order which are not of the type studied in this paper appear in the modern theory 
of pseudo-differential operators |55]. However, for the purposes of this paper, the symbolic 



calculus described in the classical papers |51] and will be sufficient, since our main 
focus is on the determination of the number of parameters on which the analytic local 
solutions of the initial value problem could possibly depend. It is a bit surprising that the 
physics community has apparently not been aware of this classical mathematical literature, 
given the simplicity of the mathematical formalism which relies only on some basic results 
from complex analysis and the theory of integral transforms. One of our primary goals in 
the current note is to bring this mathematical literature to the attention of the physics 
community and to apply the formalism to certain equations of particular physical interest. 

Since differential equations of A^-th order (in the time derivative) require N initial 
data it is sometimes reasoned that equations of the form ( |1.3| ) admit a unique solution 
only once infinitely many initial data are specified (as in our toy example ( |l.7D ). Such a 
situation would pose serious difficulties for any physical application: by suitable choice of 
infinite data one could freely specify the solution 0(t) to arbitary accuracy in any finite time 
interval At |^9| . In such a situation the initial value problem completely looses predictivity. 
Fortunately, it is not generically the case that differential equations of infinite order admit 
infinitely many initial data. We will show that for free field theories every pole of the 
propagator contributes two initial data to the solution of the field equation. This result is 
simple to understand on physical grounds since each pole of the propagator corresponds to 
a physical excitation in the theory and, on quantization, the two initial data per degree of 
freedom are promoted to annihilation/creation operators. 

In the context of quantum field theory the question of counting initial data is intimately 
related to the question of whether the theory suffers from the presence of ghosts - quantum 
states having wrong-sign kinetic term in the Lagrangian. The presence of ghosts signals a 
pathology in the underlying quantum field theory since these states either violate unitarity 



-4- 



or else carry negative kinetic energy and lead to vacuum instability |57]. One must worry 
about the presence of ghosts in theories which give rise to equations of the form ( |1.3| ) 
since the addition of finitely many higher derivative terms in the Lagrangian generically 
introduces ghost-like excitations into the theory. As we shall see later on, this is not 
necessarily the case in nonlocal field theories containing infinitely many higher derivative 
terms. A related worry is the presence of the Ostrogradski instability [^4| (see |48, 49 1 
for a review) which generically plagues finite higher derivative theories. The Ostrogradski 
instability is essentially the classical manifestation of having ghosts in the theory. Later 
on we will elucidate more carefully the relationship between these issues. 

Though it is possible to constuct nonlocal equations which evade the Ostrogradski 
instability, we will show that neither the dynamical equations of p-adic string theory (|1.1|) 



nor SFT (1.2) belong to this class. However, we illustrate a simple nonperturbative pre- 
scription for re-defining these theories in such a way as to evade such difficulties. This 
prescription restricts the theory to frequencies uj within some contour C in the complex 
plane and may be analogous to placing an ultra-violet (UV) cut-off on the theory. 

The plan of this paper is as follows. In section || we consider linear ordinary differential 
equations of infinite order, laying down the necessary formalism for counting initial con- 
ditions. In section |^ we generalize this analysis to partial differential equations of infinite 
order (in the case that the derivatives appear in the combination □ = —d^ + V^) showing 
that the initial data counting has a transparent physical interpretation and commenting 
also on the failure of the Ostrogradski construction. In section ^ we apply our prescrip- 
tion for initial data counting to some particular infinite order differential equations which 
appear frequently in the literature. In section ^ we discuss the nonlinear problem and 
illustrate the application of our formalism to nonlinear equations by studying the equation 
of p-adic string theory using a perturbative expansion about the unstable vacuum of the 
theory. We present our conclusions in section ^. Appendix A reviews the relation between 
( |1.1| ) and a certain nonlinear Fredholm equation; appendix B gives some technical details 
concerning our conventions for the Laplace transform and appendix C applies our initial 
data counting to open-closed p-adic string theory. 



2. Linear Ordinary Differential Equations of Infinite Order: Counting the 
Initial Data 

Before specializing to particular equations of immediate physical interest we first develop 
the general theory of linear ordinary differential equations (ODEs) of the form 

fidt)cPit) = Jit) (2.1) 

where the function f{s) is often called the generatrix in mathematical literature. Equations 
of the form ( ^.l]) are closely associated with both Fredholm integral equations and this 
relation is illustrated in appendix A for the case of the p-adic string. There are various 
prescriptions which one might use to define the action of the formal pseudo-differential 
operator f{dt). In the case that f{s) is an analytic function one might represent it by the 
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convergent series expansion 



so that (U^) can be written as 



n! 

n=0 



n=0 

In the case where f{s) is non-analytic however, some alternative definition is required. 
(Notice that the possibility of non-analytic generatrix is not purely academic. The zeta- 
strings model of p^ ] involves a nonanalytic kinetic function.) A natural prescription is 
to define the pseudo-differential operator f{dt) through its action on Laplace transforms.^ 
Let us assume that we are given a function (j)[t) which admits a Laplace transform (j){s) 
defined by 

m = ^^£dse''4>{s) (2.3) 

valid for t > where C is a contour to be specified (see appendix B for details^). For 
t < this expression does not, strictly, apply and (p{t) should be taken to vanish (see again 
appendix B for more details). It is natural to define f{dt)4>{t) by 

where we drop the additional term involving (0) (see appendix B) which can be absorbed 
into the arbitrary coefficients in the solutions of ( |2.lD (this will become clear shortly). If 
we make the natural choice of taking C to be an contour which encloses all the poles of the 



integrand in (|2.3| ) then (2.4) reproduces the infinite series definition (2.2) for f{s) analytic. 



Moreover, this definition reproduces the one used in |58]. This is the choice of C which is 
mathematically well- motivated and we will make this choice for the most part. However, 
later on we will consider a particular alternative choice of C which is motivated by physical 
considerations. 

We will be particularly interested in equations of the form ( |2.1| ) for which the generatrix 
may be cast in the form 

M 

f{s)=j{s)ll{s-s,Y^ (2.5) 

i=l 

with 7(5) being everywhere non-zero, which implies that /(s) has precisely M zeros at the 
points s = Si, the i-th zero being of order r^. For simplicity we assume that all the are 



positive integers, though this assumption is relaxed in subsection 2.4. We further assume 
that |si| < |s2| < • • • < \sm\- The function 7(5)"^ is nonsingular on the disk |s| < \sj\j\ but 
is otherwise arbitrary. It is useful also to introduce the resolvent generatrix f{s)~^ which 
has simple poles at the points s = Si, the i-th pole being of order r^. 

^One might instead imagine using a Laurent series definition, however, such a definition would be 
inherently ambiguous because /(s) would admit different Laurent series expansions, each valid in a different 
annulus about s — 0. 

^One could require for example that 4>{t) be entire, of exponential type. 
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2.1 The Homogeneous Equation 

We first assume that we are given a solution of in the case that J(t) = 0, and that this 
solution admits a Laplace transform (j){s) defined by ( ^.31 ). Substituting (2.3) into (p.l|), 
yields 

±,£dse''f{s)4>is) = (2.6) 

What is the most general function (p{s) which satisfies this condition? From the Cauchy 
integral theorem we know that this condition will be satisfied only if the integrand is 
analytic everywhere in a neighborhood of the interior of C so that (p{s) may have simple 
poles at the points s = Si, the i-th pole being of order rj or less. Thus the most general 
solution can be written in the form 



4>^s) = ^yy , ' (2.7) 

The factor of 7(5) is included for convenience and an additive analytic function (which 
would not alter the configuration-space solution (/>(t)) has been omitted.^ The solution 
( |2.7| ) contains a total of arbitrary coefficients A^j^ where 

M 



N = J2rt (2.8) 



i=l 



In other words, N counts the zeros of f{s) according to their multiplicity. The N free 

(i) 

coefficients Aj will ultimately serve to fix N initial conditions for the initial value problem 



corresponding to ( |2.lD . We now insert the solution ( p.7| ) into ( |2.3[ ) and perform the ds 
integration (with C a contour which encloses all of the points {si}) using the Cauchy 
integral formula 

2m Jc {s - Si)! (j - 1)! 

(valid for h{s) analytic inside C and j > 0). The resulting solution (j){t) in configuration 
space takes the form 

M 



i=l 



where each Pi{t) is a polynomial of order rj — 1 

j=i 

The coefficients {p^j^ } are arbitrary (refiecting the fact that A^j^ were arbitrary) and will 
serve to fix N initial conditions ^'•"'''(O) for n = 0, • • • , — 1. 



^The additional terms excluded from (^.4|) (see appendix B) can be absorbed into the definition of 
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2.2 The Particular Solution of the Inhomogeneous Equation 

Having determined the solution of the homogeneous equation (|2.1|) we now consider the 



situation where J{t) / and focus on the particular solution due to the source J{t). We 
assume that the source term has a Laplace transform J(s) defined by 



J{t) = ^ (j^ dse'^J{s) (2.11) 



27ri 



and we continue to employ the Laplace transform (p{s) (|2.3| ). Then equation {2A) takes 
the form 



1 
27ri 



ds 

C 



fis)^{s) - J{s)\ = (2.12) 
so that the particular solution is cj}{s) = J{s)/ f{s) or, in configuration space 



From (|t|) it is clear that the resolvent generatrix is the Laplace-space Green function 

G{s) = f{s)-\ 

2.3 The Initial Value Problem 



It is natural to expect that in order to obtain the most general solution of (2.1) we should 



append to ( 2.13| ) the solution of the homogeneous equation (2^). One may prove the 



following theorem, which we now state without proof but whose verisimilitude should be 
clear from the preceding analysis (for a detailed and elegant proof, see [|53|). 

Theorem: Consider the linear differential equation of infinite order ( \2. 1\ ) with the 
function f{s) analytic within the region \s\ < q where q is a given positive constant or zero 
and further suppose that f{s) may be written in the form ^2^. If J{t) is a function of 
exponential type not exceeding q then the most general solution (p{t), subject to the condition 
that it shall be a function of exponential type not exceeding q, may be written in the form 

where Pi{t) is an arbitrary polynomial of order r j — 1 and J{s) is the Laplace-space source. 
The first term in ( 2.1^ ) corresponds to the particular solution due to the source J{t) and 



the second term is the solution of the homogeneous equation. If the generatrix f{s) does 
not vanish inside \s\ < q then the latter solution is identically zero. 

Armed with the solution ( |2.14| ) the interpretation of the initial value problem for 
(|2.1|) is clear: the solution contains N arbitrary coefficients pf^ which serve to fix 



initial conditions (f>^"'\0) for n = 0, • • • ,N — 1. This result has a transparent physical 
interpretation in terms of poles of the propagator, which we shall return to in the next 
section. It is important to note that the significance of this theorem is that ( p. 141) provides 
the most general solution of (|2.1| ) . Returning to the motivational example (|1.7] ) we see that 
( |1.8| ) is, indeed, the full solution because the zeros of the Airy function are order unity. 
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An interesting consequence of this theorem, whose significance we shall return to later 
on, is that for J{t) = the solution ( |2.14| ) is identical to the solution of the equation 
f{dt)(p{t) = where 

m = 44 = - '^y- 

7(s) l\ 

In other words, for the homogeneous equation the dynamics are completely insensitive to 
the choice of 7(5) and the solutions are completely determined by the pole structure of the 
resolvent generatrix. For finite M this implies that the dynamics of the full infinite order 
differential equation (provided it is linear and source-free) will be identical to some finite 
order differential equation. 

It is worth commenting on the generality of our results. From the perspective of 
mathematical rigour the only serious caveat of our analysis is the assumption that the 
sources J{t) and solutions (p{t) are analytic, with the asymptotic behavior needed for their 
Laplace transform to exist (which is not, in general, guaranteed). Nonanalytic solutions 
have been discussed for p-adic string theory in |^0| though it is not entirely clear how to 
interpret these solutions physically. Naively, one might expect some components of the 
stress tensor associated with such solutions to blow up at some finite time. However, Q 
showed that a particular nonanalytic solution in SFT could sensibly be regularized to yield 
a well-behaved stress tensor. 



2.4 Equations Involving Fractional Operators 



In writing ( |2.5| ) we have assumed that the zeroes of the generatrix are of integer order (we 
have assumed that all of the rj are integer) and we have excluded the case of fractional dif- 
ferentiation (/(s) = with Q non-integer). In this subsection we relax those assumptions. 
Though such equations are less common in applications to string theory and cosmology, 
they constitute a large class of nonlocal equations which admit well-posed initial value 
problem with finitely many initial conditions and may be of some interest for this reason. 
The reader more interested in the equations {lA, |1.2| ) may wish to skip to the next section. 



2.4.1 Roots of Non-Integer Order 

To illustrate the effect of taking non-integer in (^]B) let us consider the simple equation 



+ m2 ^{t) = 



(2.15) 



If we seek solutions (j){t) which admit a Laplace transform (j){s), we obtain 

1 



2iTi 



ds 



C 



\/s2 _^ ^2 ^[s)e'^ 







which implies that 



A 



(2.16) 



(2.17) 



-9- 



up to an additive analytic function which does not alter the configuration-space solution. 
Taking the inverse Laplace transform using the formula 



£-1 



1 



Jo(mt) 



(where J^, is the Bessel function of the first kind) we have 

4){t) = AJoimt) 



(2.18) 



(2.19) 



So that equation ( |2.15| ) admits only one initial condition. 
It is also straightforward to consider 



which implies that 



The general solution is 



1 

27ri 



C 



{d^ + (l>{t) = 



ds 



+ 



(s2+ ^2)1/2 (^2+ ^2)3/2 



(2.20) 
(2.21) 

(2.22) 



with o, b arbitrary. We must now perform the inverse Laplace transform. The first term is 
trivial using ( p. 18 ). By differentiating ( 2.18| ) with respect to m we find that 



(s^ -|- m^)^^^ 



m 



-Ji{mt) 



so that the configuation-space solution of ( |2.20| ) is 

4){t) = AJo{mt) + B{mt)Ji{mt) 

which admits two initial conditions. 

It is straightforward to handle higher order equations of the form 

(52 + m2)"/'0(t) = O 

with n an odd integer. The Laplace-space solution is 



(2.23) 



(2.24) 



(n-l)/2 



i=0 



(s2 _^^2)(2i+l)/2 



which contains a total of (n -|- l)/2 arbitrary coefficients. The inverse Laplace transforms 
may be performed by successively differentiating ( p. 18 ) with respect to m and using the 
well-known Bessel function identity 



J'Ax) = -Ju+i{x) + -Ju{x) 

X 
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In the tachyonic case m? = — /^^ < one simply replaces the Bessel functions Jj/ with 
modified Bessel functions /jy to obtain real-valued solutions. We find that the equation 

^IW^mt) = (2.25) 

has solution 

(^(t) = AIo{^lt) (2.26) 

while 

{df-f,Y'<Pit) = (2.27) 

has solution 

(Pit) = Aloifit) + B{fit)Ii{fit) (2.28) 
2.4.2 Fractional Differentiation 

Another interesting case occurs when the generatrix involves noninteger powers of the 
derivatives, such as /(s) = with a some arbitrary complex number. In this case clearly 
the definition ( p.2| ) does not apply. The attempt to assign sensible meaning to the symbol 
for noninteger a has a long history in mathematical analysis and there are a number 
of definitions in the literature |5^, For our purposes, the most useful definition seems 
to be the Liouville fractional derivative. To motivate this definition we first consider the 
problem of defining fractional integration. We define ordinary integration by 

(/</>) (t) = Tdtxt') 







{I'm = / dt'imt') 





(The choice of lower limit of integration is purely conventional. Different choices will lead 
to different definitions of fractional differentiation.) The Cauchy formula for repeated 
integration gives 

(/»(t) = ^^^^ tdt'it - t'r-'mdt' (2.29) 







A natural continuation to noninteger n is 

1 



(r </.)(*) = / * dt'it - t'r~'mdt' (2.30) 

a — 1 ! 







valid for t > 0, Re(a) > and where the symbol (a — 1)! should now be interpreted as 
the factorial function.^ This definition is commutative and additive I"(I^(/)) = /^(I"0) = 
It is now natural to define fractional differentiation by differentiating ( p. 30]) 

-J (/"-»(*) 

1 ( dY f\, Mt') , , 

d^i ^1 2.31 

/ 4. 4./\rv— rj-H \ ' 



n-a-\)\\dt) Jo (t-t')°""+^ 



*That is to say [a — 1)! = r(Qf) (where F is the gamma-function) for non-integer a. We do not use 
the gamma-function notation explicitly to avoid confusion with another function which we denote by r{z) 
which will be introduced in the next section. 
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Of course for integer a equation ( 2.31 ) reproduces the usual derivative. This definition also 
has the following nice properties 15^ 



{D'^I'^(^){t) = (l){t) if a > 
[D^I'^(^){t) = (/°-^(/>)(t) if a > /3 > 



(/?-a- 1)! 
The definition ( |2.30| ) acts in Laplace space as 



if Re(a) > 0, Re(/3) > 



S 0[S) 



where I is a natural number such that I — I < a < I and the constants {di} are 



(2.32) 



(2.33) 



The result (2.32) is exactly of the form d2 



2.4.3 Fractional DifFerential Equations 

There is a large literature on solving equations involving the operator (2.31) which is on 
rather firm mathematical footing (both linear and nonlinear equations have been studied). 
For illustration we consider only the simple prototype equation 

(L>°0)(t) - m"0(t) = J{t) (2.34) 

We first consider the homogeneous equation J = 0. Taking the Laplace transform using 
(2.32) and solving for <p{s) we have 



Now, inverting the Laplace transform we find the solution 



where 



'jit) = 



(2.35) 



where E^^is is the generalized Mittag-Leffler function (see, for example, [§9|). The funda- 
mental solutions {(j)j} satisfy 

(I)°->,)(0) = for i,j = l,--- ,l;i>j 



(i?°-V,)(0) = 1 for i = l, 
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It can be verified that in fact any summation 



(l){t) = Y,aj(pj{t) (2.36) 



witli / arbitrary coefficients {aj} yields a solution to ( p.34 ). We see, then, that equation 



(2.34) admits a well-posed initial value problem with / initial conditions. In particular, for 



1 < a < 2 equation ( p. 34 ) admits two initial conditions. 



It is also straightforward to consider the inhomogeneous equation, following the ap- 



proach of (2.13). In fact, one may define a fractional Green function 



1 f e"^ 

GJt) = i ds (2.37) 

^' 2-711 Jc s'^-m°' ^ ^ 



so that the particular solution is 



0(i) 



/ dt'Ga{t-t')J{t') (2.38) 

^0 



which can be proved to provide a unique solution to (2.34) with (j){0) = 



2.4.4 Properties of Fractional Operators 

It is worth commenting on the fact that our definition of the operator {df + m^)^/^ and 
also the Liouville^ definition of the fractional derivative -D" lead to some a-priori unex- 
pected properties, the most notable of which is that the degrees of the fractional derivative 
operators which we have defined are not, in general additive. In other words, D"' does 
not satisfy D°^D^ = for general (non-integer) a, j3 |^^. As simple illustration of 

this fact is given by considering the function (f){t) = which satisfies D^f^cf) = 0, while 

D^cf) = dtcj) = —1/2 7^ 0. A similar property can be inferred in the case of the operator 
{df -|- m^)^/^ by noting that ( |2.15| ) does not have solutions cos(mt), sin(mt). Indeed, it is 



straightforward to see that the solution ( 2.19 ) of ( 2.15 ) is not a solution of the equation 
((9| + m2)0(t) = 0. 

Given that, using our conventions, acting twice with the operator (9^ -|- m?Y^'^ does 
not necessarily return the same result as acting once with the operator {di + m^), the 
reader may wonder whether the definition that we have employed is the most sensible one. 
We would like to argue that, although the definition we have adopted has the unpleasant 
property that the degrees are in general not additive, it is the most physically reasonable 
approach. To see this, we consider an alternative definition of (9f +m^)^/^ and contrast this 
with ours. As an alternative to our approach let us consider defining f{dt) = \/&f+rn?' 
by an infinite series of the form 



+ ^2 = 5t + Qo + Q-i + Q-2 + Q-2 + ■■■ (2.39) 



''Our comments actually apply for nearly all definitions of fractional differentiation which are considered 
in the mathematics literature. See ISyl for further details. 
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Substituting this into the defining identity 



and applying the formal identity 

oo 

to solve for the coefficients Qo 



=0 



Q-ii Q-2, Q-3 ■ 



d? + m'^ = dt + 



m 



we arrive at 
1 



■ m 



87' + ... 



(2.40) 



(2.41) 



(2.42) 



It is ensured by construction that acting twice one some test function with this operator 
will return the same result as acting once with df + . 

In order to see how the definition (2.42) differs from our Laplace-space definition (| 
we consider acting with ( p. 42 ) on a function of the form ( |2.3| ). The result is 



ds e 



St 



C 



m 

S + —S- 



1 4 - 

— ms 



+ 



(2.43) 



The reader will recognize that the series in the square braces coincides with the Taylor 
series for + m? about m = 0. In the region \m/s\ > 1 (the high energy part of the 
phase space) the series converges to + m?. However, in the region \m/s\ < 1 the 
infinite series fails to converge. It is due to the behaviour of the infinite series in this low 
energy part of phase space that ( 2.43| ) and (2^) do not yield the same results. (Notice that 
the series fails to converge at the zeroes of the generatrix = —m? which is not surprising 
because these are branch points of the function Vs^ + m'^-) 

The advantage of the definition ( 2.42 ) over our approach (2^) is obvious: the series 
definition preserves the property ( 2.401 ) which one expects for the square root of an operator. 
However, we will argue that our approach seems more sensible from the perspective of 
constructing quantum field theories. Since the series ( 2.42| ) fails to converge at \s\ <C \m\ 
this implies that there is no well-defined low energy limit. In fact, using this definition, 
it does not seem possible to recover the usual Klein-Gordon Lagrangian at low energies. 
Moreover, the computation of scattering amplitudes will involve integrating the propagator 
over all momenta. Using the definition ( 2.42| ) the low energy part of the phase space 
integrals will necessarily be ill-defined. Finally, we point out a more heuristic problem 
with the series definition (2.42). This definition fails to converge at the curical point 
= — and thus it seems necessary to adandon the usual interpretation of the poles of 
the propagator as physical states in the associated quantum field theory. 

We should stress that our definition ( |2.4| ) is just that and that it may be interesting 
to consider other definitions such as ( 2.42| ) which do not lead to the unattractive property 
which we have discussed above. We would also like to stress that that precisely the same 
unattractive feature arises for all of the definitions of fractional differentiation Z?" which 
are studied in the mathematics literature. Thus, it is clear that one need not view this 
complication as prohibitive to any potential physical application of such operators. The 
reader who disagrees with this statement may wish to view the analysis of this subsection 
as a motivation to avoid nonlocal theories involving fractional operators. 
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3. Physical Interpretation of the Result 



3.1 Linear Partial Differential Equations of Infinite Order 

The analysis of section |2| is actually somewhat more general than what is required in order 
to study ( |1.1| , p^ ). For differential equations of infinite order which arise from Lorentz 
invariant field theories the time derivatives dt must always appear within the d'Alembertian 
operator □ = —d^ + V^. Hence we would now like to apply the preceding analysis to linear 
partial differential equations (PDEs) of infinite order of the form 

F(n)0(t,x) = J(t,x) (3.1) 



We refer to the function F{z) in ( p.lD as the kinetic operator, which is closely related to 
the generatrix. We will be particularly interested in kinetic operators which can be written 
in the form 



TV 



F{z) = T{z)J{{-z + m^,) (3.2) 



i=l 



where V{z) ^ contains no poles in the complex plane, but is otherwise arbitrary.^ With the 



ansatz (3.2) equation ( p.l| ) describes N physical states with masses {rm}. For simplicity 
we assume the rrii to be non-degenerate, however, it is simple to drop this restriction. 
(Including degenerate masses corresponds to choosing some of the in equation (2.5) to 
be different from unity.) 

We now proceed to construct the particular solution of (|3.lj ). Assuming that the field 
(j){t, x) can be expanded into Fourier modes as 

equation (|3.1| ) becomes 

F{-dj - k^)Ut) = Mt) (3.4) 

where fc^ = k • k, /c = \/A? and 

Mt)= [ -^^e-^'^-V(t,x) 



(27r)3/2 

Equation ( ^.4| ) is of the form ( |2.lD where the generatrix is 

fis) = Fi-s'' - k") 

N 

= T{-s' - k') n + iiJf^) {s - i4^) (3.5) 



i=l 



SO that the resolvent generatrix has two poles for each pole of the propagator. In ( ^.5]) we 
have defined 



®r(2) should not be confused with the well-known gamma- function. 
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Following we see that the particular solution of ( p.4| ) may be written as 

Ms) 



St 



F( 



so that Laplace-space Green function is Gk{s) 
propagator is 

= Gkiiio) 



2\-l 



1 



F( 



(3.7) 

The momentum-space 
(3.8) 



with p'^ = —u? -I- . 

We now proceed to solve (^) for the homogeneous case J = 0. Since the resolvent 
generatrix has 2N poles (of order one) we expect the solutions to contain free coefficients 
for each fc-mode (two for each physical degree of freedom). It will be convenient to write 
these 2A^ free coefficients as the real and imaginary parts of N complex numbers a^*\ 



In solving (3J.) we wish to apply an additional constraint on the solutions which was 
not implied by the analysis of section (^). Namely, we demand that the solutions (/>(t,x) 
be real valued. The condition 6 = 6* translates into the constraint 



(3.9) 

on the Fourier modes. The general solution of (|3.l| ) consistent with the reality condition 
(|3.9|) then must take the form 



N 



ek(t) = E^k^w 



it) 



1=1 



(t)+a(_t</'l^(t) 



(3.10) 



(3.11) 



Each is the solution corresponding to the i-th pole of the propagator (with mass m'f) 
and </>k(i) are the mode functions (which have been constructed in such a way that (pi^ and 
(p^ are linearly independent). The configuration-space solution can be written as 



(/)(t,x) 



N 

E 

i=l 



(27r)3/2 



a!^4\t)^''-^ + c.c. 



(3.12) 



where "c.c. 



denotes the complex conjugate of the preceding term. In the classical solu- 
tion the coefficients , a^''* serve to fix 2N initial data 0(0, x), 0(0, x), • • • , ^^^a^-^) (0, x). 
However, in the quantum theory these coefficients are promoted to annihilation/creation 



operators. It is clear that in this context the solution ( 3.12 ) describes physical degrees 
of freedom, one for each pole of the propagator. 

The solutions evolve differently depending on the value of mj. There are three 
qualitatively distinct cases: 

1. Stable modes: mf > 0. This case corresponds to real u!^\ The solution of ( |3.4| ) takes 
the form (3.10, 3.11) with mode function 



■ (i). 



2uj 



(3.13) 
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(The normahzation of the modes t/)^*^ is purely conventional since the constant pre- 
factor may be absorbed into the arbitrary coefficients ai^ ■) 



2. Tachyonic modes: mf 



< 0. In this case io^^ is real for k'^ > 



^ and pure 



imaginary for k^ < fif. We consider only the latter modes (the instability band) 
since the former possibility is identical to the previous case. For tachyonic modes 
within the instability band the solution of (3^) takes the form ( 3.10| , 3.11 ) with 
mode functions 



1 



2\ 2n 



(0 



(3.14) 



where ri 



? — k'^ is real-valued. 



3. Poles with complex mass. Taking m~ to be some arbitrary complex numbers having 
nonvanishing imaginary part the frequency can be written as 



a 



A crucial point is that for a CPT invariant theory the complex-mass states must 



always arise in conjugate pairs. Hence we restrict ourselves to the case where lv 



'^k^ — iP^!' is also a pole. Real-valued solutions can be obtained by appropriately 
superposing the particular solutions corresponding to the poles and ivt as 



+ 



2\uj 



2\uj 



-la, t 



(3.15) 



In general the summation (3.12) will contain all three types of modes ( p. 13 , 3.14 
3.15| ), although any particular class could be consistently projected out by altering the 
prescription for drawing the contour C (see the discussion below equation ( |2.4D ). For 
example, one could exclude the complex mass solutions by taking C in (|2.3|) to encircle 



only poles of the resolvent generatrix which lie on either the real axis or the imaginary 
axis in the complex plane. Later on, we will employ this particular choice of C to render 
both p-adic string theory and SFT ghost-free. Though such a prescription is well defined 
mathematically, it is not clear how to interpret the resulting truncated theory from a 
physical perspective. However, one might imagine taking the contour C as a part of the 
definition of field theory so that different choices of contour yield different theories with 
different mass spectra. Since this prescription restricts the theory to complex frequencies 
u! inside C it is, in some sense, analogous to putting a UV cut-off on the underlying 
field theory. This procedure of projecting out unwanted states by suitable choice of C 
seems reasonable at the level of effective field theory (and hence should be acceptable for 
applications to cosmology), however, it is less clear if this is sensible in string theory which 
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is supposed to be UV complete. On the other hand, notice that the particular contour 
which we have discussed still allows for infinitely large frequencies. Rather, what we are 
constraining is the direction in the complex plane in which the frequency can be made 
large. 

3.2 Generalization to Curved Backgrounds 

Though our analysis has been restricted to 3+1-dimensional flat Minkowski space, it should 
be clear that these conclusions readily generalize to D-dimensions. In curved backgrounds 
one might consider generalizing our analysis to equations of the form 

F(n3)0(t,x) = J(t,x) (3.16) 

where Og = g^'^'V^'Vi, is the covariant d'Alembertian operator and we take F{z) of the 
form ( |3.2| ). In the homogeneous case J = a solution may be constructed by taking 

N 

(P{t,x) = ^Mt,^) (3.17) 

1=1 

where each (j)i{t,x) obeys an eigenvalue equation 



g<pi = mi(t>i (3.18) 

(this approach of taking eigenfunctions of the d'Alembertian was employed to study the 
cosmology of nonlocal theories in |32] and also in [27, |3^, ^). Each solution of ( 3.1^ ) 



(assuming that solutions exist) should admit two initial data and hence the full solution 
( |3.17 ) admits 2A^ initial data, as in our previous flat-space analysis. However, it is not so 



clear if this analysis can be generalized to include inhomogeneous equations of the form 
(3.16). Any specific analysis will require detailed knowledge of the pole structure of the 



propagator which is, in general backgrounds, a highly nontrivial task. 

It is, however, straightforward to generalize our discussion to include homogeneous 



solutions in de Sitter space. For a kinetic function of the form (3^) the de Sitter space 
generatrix is 

/(s) = r(-s2 - 3Hs) Yl (s^ + SHs + ml) (3.19) 

n 

(the restriction to homogeneous solutions means that we are considering only k = 0). There 
are two zeros for each mass m„ 



(±) _ _3H 



4m2 



(3.20) 



It is interesting to consider the limit of large Hubble friction S> in which case the 
solution takes the form 

(t>{t) ^ 0oe-^^* + a„e-™"*/(^^) (3.21) 

n 

In the limit i7 — > oo the first term damps to zero quickly while the second term approaches 
a constant. The limit of large Hubble friction is relevant for studies of inflation in nonlocal 
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theories since during inflation m? /H"^ ~ C(^) where m is the characteristic mass scale in 
the problem and r/ <C 1 is a slow roll parameter. In this case the almost-constant term in 
(3.21) coincides with the slowly rolling solution. 

3.3 Partial Differential Equations Involving Fractional Operators 



Using the results of subsection 2.4.1 it is straightforward to consider equations of the form 

(-□+m2)"/2 </.(t,x) =0 

with n an odd integer. Assuming that (p(t,x.) can be represented as a Fourier integral, we 
have, in Fourier space. 



(92 + fc2 + m2)"/'ek(t) = 



(3.22) 



Let us first consider n = 1 which has only a single mode function. For stable modes > 
we have the solution 



Ck(i) = (ok + a-k) 



1 



JoiuJkt) 



(3.23) 



where oj^ = Vk'^ + m'^- (This solution is of the form ( p.ll|) though it only contains one free 
coefficient because 0k and (;/)_k are not, in this case, linearly independent.) At late times 
LOkt ^ 1 the solution (3.23) undergoes damped oscillations 



Ck(i) = (flk + a\] 



For tachyonic modes = —fi^ < we have the solution 



Ck(i) = (ok + a^k) 



1 



/o(Ofct) 



(3.24) 



where = \/ — k'^ and we are assuming that k < fi. As one would expect the solution 
(3.24) grows exponentially 



cut) = (ak + a-k) 



A 



2Qi.V7rt 



at late times Qkt ^ 1- The case m = has been studied in 2 + 1-dimensions 



where 



it has applications to bosonization [61|. (Equations involving fractional powers of the 
d'Alembertian have also been studied in |62].) This special case admits a canonical quan- 
tization and both causality and Huygens' principle are both respected ||60t] . 

We now consider ( p.22[ ) with n = 3. In the stable case m? > the solutions can be 
cast in the form ( 3.11j ) with mode function 



while, for the tachyonic case m? 



[JoiuJkt) + i{LJkt)Jl{LOkt)] 



2ujk 

we would have 

[io{nkt) + i{nkt)ii{Qkt)] 



1 



(3.25) 



(3.26) 
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It is also straightforward to consider partial differential equations involving fractional 
derivatives. Sacrificing Lorentz invariance we might consider the equation 



Z?"(/.(t, x) = \^did'(t){t, x) (3.27) 

which interpolates between the wave equation (when a = 2) and the diffusion equation 
(when Q = 1). After performing a Fourier transform with respect to the spatial variables 
( |3.27 ) becomes an equation of the form ( 2.34 ). Thus, for 1 < a < 2 equation ( |3.27| ) admits 



a well-posed initial value problem with two initial data and for < a < 1 it admits only 
one piece of initial data. It is straightforward also to add a mass term to ( |3.27| ). Equations 
of this form (and also many more general fractional PDEs, both linear and nonlinear) have 
been studied in great detail in the mathematics literature pOf . 



3.4 Ghosts and the Ostrogradski InstabiUty 

Though essentially all nonlocal theories containing finitely many higher derivatives have 
ghosts, this is not generically true of theories with infinitely many derivatives. For example, 
in the case where the propagator contains no poles it is clear from ( p. 12 ) that the underlying 



field theory has no physical excitations at all, ghost or otherwise. The question of whether 



or not ghosts are present in an infinite derivative theory was considered in |22| where it 
was shown that the theory will be ghost-free as long as the propagator contains at most a 
single pole (of order unity). However, if the propagator contains two or more poles then the 
theory will almost always contain ghost-like excitations.^ Physically this result is easy to 
understand: theories with only a single pole in the propagator describe only one physical 
degree of freedom and hence the nonlocal structure does not introduce spurious new ghost 
states. 

If we restrict ourselves to real-valued solutions of linear differential equations which 
arise from Lorentz- invariant, ghost-free theories with analytic generatrix then we are limited 
to equations of the form 

r(n)(-n + m2)</.(t,x) = J(t,x) (3.28) 

with r(z) analytic everywhere in the complex plane. In the case J = 0, as we have previ- 
ously discussed in subsection |2.3| , the solutions of this equation are completely insensitive 
to the choice of r(z). In particular, the dynamics for J = are identical to the solutions 
of the local wave equation 

(-□ +m2)(/)(t,x) = 

so that in the linear, source-free theory the nonlocal structure has no effect on the dynamics 
(with the possible exception of re-scaling the quantity which is naively identified as the mass 



of the particle). This result was also observed in |32], though it was not rigorously justified 
(this result was also observed in [^^). Equations of the form (3.25) fall into the class of 
theories referred to as trivially nonlocal in |^9[ where it was argued that the nonlocality 
can be removed by a field redefinition. (This is, of course, only true at the linear level. 



^Under very special circumstances it is in fact possible to construct ghost-free multi-pole theories 
However, we restrict ourselve to single-pole theories for the present analysis. 
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Inclusion of a nonlinear term in ( 3.2^ ) will spoil the triviality. It is, unfortunately, not 



clear whether such an addition will generically also spoil the well-posedness of the initial 
value problem.) 

One can also construct ghost-free linear differential equations by considering kinetic 
functions F{z) with a single zero and also a single pole. Such equations can be cast in the 
form ( p. 28 ) with T{z) having single pole at some point z = jj,'^ and may admit well-posed 



initial value problem, even at the nonlinear level. For example, the equation 

f{dt)m = um)] (3.29) 

with generatrix 

2^1 / s2x 



f{s) = — (^1 - -J (3.30) 

can be seen to be equivalent to the local equation 

{df - A^) U [^{t)] + 2AX^{t) = (3.31) 

by acting on both sides of ( p. 291 ) with the operator f{dt)~^. One expects the initial value 
problem for ( |3.31 ) to be well-posed with only two initial data. (Similar conclusions apply 



also to more general equations which involve f{dt)(t>{t) and its first two derivatives.) The 
well-posedness of (|3.29| ) is not surprising since this equation can be obtained from the 



local theory ( 3.31 ) by acting with an inverse differential operator f{dt)^^- This kind of 
nonlocality is referred to as derived in |49|] . 

Finally, we note that it seems to be possible to construct a wide array of ghost-free 



theories using fractional operators, as in subsection 3^. In the case of equation (|3.271) with 



< a < 2 we have an example of a nonlocal theory with well-defined initial value problem 
which is not derivable from some local theory, though the price is manifestly broken Lorentz 
invariance. Are equations of the form ( |3.22 ) derivable from local theories? For an equation 



of the form (— □ -|- m?Y~°'(j) = one might be tempted to define = (— □ -|- m?')^°'x so 
that X obeys the local equations (— □ -|- rTn?)x = 0. However, it is argued in [^] that such 
a procedure is ill-defined. Indeed, the transformation (p ^ x has zeroes in the complex 
plane and would not usually be allowed in a field redefinition. In the case a = 1/2 it is 



clear from ( 3.23 ) that such a procedure may count incorrectly the number of initial data. 

We have seen how the question of counting initial data is intimately related to the 
question of whether the underlying field theory contains ghosts. We now consider a re- 
lated worry: the Ostrogradski instability [44|. A cartoon of the Ostrogradski construction 
follows (for a modern review see [|^, ^^). Consider a higher-derivative Lagrangian which 
depends nondegenerately on the field and its first derivatives (so that the equation 
of motion is of order 2N in time derivatives) the Hamiltonian depends on 2A canonical 
coordinates corresponding to the 2A pieces of initial data which are necessary to spec- 
ify the solutions of the Euler-Lagrange equation. Ostrogradski's theorem states that the 
Hamiltonian always depends linearly A — 1 of the conjugate momenta. It follows that 
the Hamiltonian is generically unbounded from below and hence it is necessarily unstable 
over half the phase space for large A. In light of the previous analysis it is easy to see 
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why Ostrogradski's construction may fail for theories containing infinitely many deriva- 
tives. As long as the propagator only has one pole, then only two initial data are necessary 
to specify the solutions of the equation of motion corresponding to only two independent 
canonical coordinates, rather than infinitely many as one would conclude by naively taking 
the N ^ oo limit of Ostrogradski's result. For nonlocal theories of infinite order one should 
instead construct the Hamiltonian using the formalism of [45, 46] which allows for a more 



transparent identification of the physical phase space of an infinite derivative theory and 
reduces to the Ostrogradski construction in the finite derivative case. 

It is interesting to note that for theories with more than one pole in the propagator 
(theories which have ghosts) turning on large Hubble friction seems to ameliorate the 
problem somewhat by slowing the evolution of the spurious modes (see equation (|3.21|) ). 
Of course, Hubble friction does not actually solve the problem by rendering the Hamiltonian 
bounded from below, it merely slows down the higher derivative instability. Obviously in 
the case of an infinite spectrum of masses {mf^} where |m^| — > oo as n — > oo then one may 
always find some level so that > for finite H. 

The failure of the Ostrogradski construction for infinite derivative theories has previ- 
ously been noted in the math literature, though not using this language. In (see also 



1 51]) it was noted that the naive procedure for writing the A^-th order differential equation 



dt ' dt"^' ' dt^ J 

as a system of A^ first order equations 

= fn it,4>,(pi,4>2, ■ ■ ■ ,(t>N) , n = l,2,--- ,N 

dt 



fails when N = oo. Though there exists no general no-go theorem which excludes the 
possibility that there exists some alternate procedure for writing an infinite order ODE as 
an infinite system of first order ODEs^, it is very unlikely that such a procedure can be 
found. 

It is interesting to note that one would obtain quite different conclusions from what we 
have discussed for the solutions of ( 3.2§| ) if one were to truncate the full kinetic function 



F{z) at some large but finite order in derivatives, A^. Consider the truncated kinetic 
function 

F(")fO) 

Ftrunc(^) = V (3.32) 

^-^ n! 

n=0 

for A^ » 1. S ince -^trunc(-^) is ^-n A^-th order polynomial in z it will inevitably contain 
some large number of spurious zeros which were not present in the full kinetic function 
F{z). Because the propagator in the truncated theory contains a large number of spurious 



*Such a situation would be reminiscent of multi-Hamiltonian systems where the second order equation(s) 
can be written as a first order system in different ways, depending on the choice of symplectic structure 
and Hamiltonian. 
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poles it follows that the solution (j){t,x) will contain a large number of modes which are 
not present in the true theory ( p.28[ ). Inevitably some of these modes will correspond to 
ghost-like degrees of freedom so the solution will be unstable. Evidently, infinite derivative 
theories only make sense when the full nonlocal structure of the theory is included. 



4. Some Linear Equations of Physical Interest 

We now apply the formalism developed in the previous section to some particular prob- 
lems of physical interest. Though we do not uncover any new solutions, our analysis places 
previous literature applying nonlocal field theories to cosmology [32, pi, p3] in the math- 



ematical context developed in this paper, shedding light on the nature of the initial value 
problem for these equations. 

We proceed by first linearizing the non-linear field equations under consideration 
around known solutions, and then solving the linearized equations using our formal op- 
erator calculus. While we don't derive the error estimates that would make our use of the 
linearized equations fully mathematically rigorous, we believe that this step is justified in 
the context of the applications of these field equations to cosmology. In the case of nonlocal 



hill-top inflationary models |32, 34 1 the background dynamics which are interesting for 
inflation occur very close to the unstable maximum of the potential and the COBE normal- 
ization ensures that inhomogeneities are very small. Thus it is quite natural to expect that 
one should obtain accurate results by linearizing the equations of motion about the false 
vacuum. Further support for this approach of perturbing about some constant solution 
is provided by the fact that our linearized solution of the p-adic string equation (below) 
essentially reproduces the leading term (at early times) in the expansion in exponentials 
used in As we shall see in section ^, working to higher order in perturbation theory 



will reproduce the full solution of |13]. Since the solutions of [13| were verified using a 



number of non-trivial consistency checks, we can be fairly confident that errors are small. 
4.1 p-adic Strings Near the False Vacuum 

The equation of motion for the tachyon field in p-adic string theory is ( |1.1| ) . This equation 
has the constant solutions (p = 0,ibl for odd p, and (j) = 0,+l in the case of even p, 
corresponding to critical points on the potential. The solution (p = 1 (and <J) = —1 for odd 
p) represents the unstable tachyonic maximum of the potential while (j) = represents the 
stable (true) vacuum of the theory. 



We first consider equation (1.1) linearized about the false vacuum which physically 
corresponds to studying small tachyonic fiuctuations which encode the instability of the 
D25-brane in p-adic string theory. Taking 

(/)(t, x) = 1 + 5(/)(t, x) (4.1) 

and linearizing (^]^) in Sep we find 

p-^/"^ -p\5(p{t,:ic) =0 (4.2) 
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The kinetic function 

F{z) =p''/^ -p (4.3) 



has infinitely many zeros at z = z„ where 



-2 + - 4.4 

mp 



for n = 0, ±1, ±2, • • • (equivalently there are infinitely many poles of the propagator at 
these points). Thus the solution of (|4.1| ) admits infintely many initial data. The solution 
(50 contains a single tachyonic mode with mass-squared —2 in string units. The behaviour 
of this mode is described by ( p. 14 ) with = 2. This is the solution which would be 



expected on physical grounds. However, there are also an infinity of complex-mass solutions 
of the form ( p.l5| ). These extra modes have equally spaced mass-squared and, perhaps, 
correspond physically to the infinite tower of massive states of bosonic string theory. It 
is not clear how reasonable this interpretation is since the (perturbative) Hamiltonian is 
unbounded from below |46].^ One could render the theory ( |l.lD physically sensible by 



simply projecting out the unwanted complex-mass poles by taking the contour C in 
to enclose only the poles along the real and imaginary axes in complex plane. We have 



discussed this prescription previously, at the end of subsection 3T. The fact that (|LJ) 
can be rendered stable^^ by suitable choice of C is consistent with [^] where it was shown 
that the Hamiltonian becomes bounded from below once suitable asymptotic conditions 
are imposed on the tachyon field. 

4.2 p-adic Strings Near the True Vacuum 

We now consider equation (|1.1| ) linearized about the true vacuum (/> = 0. Writing 

x) = + x) (4.5) 

and linearizing equation (IT) in 6(p (assuming p ^ 1) we have 

p~°/^6(t>{t,x) =0 (4.6) 

The kinetic function 

F{z) = p-'/^ (4.7) 

has no zeros in the complex plane. Thus equation (^j^) has no nontrivial solution. This 
is to be expected on physical grounds since the tachyon vacuum should contain no open 
string excitations. 

This analysis does not rule out the existence of nonperturbative solutions in the true 
vacuum of the theory. Indeed, [p!^ ] found anharmonic oscillations about (f> = using numer- 
ical methods. The anharmonic oscillations of |13| appear to contain two free parameters 
(the frequency to and a phase shift which sets the origin of time) . 



^Recall, however, that we are perturbing about an unstable maximum of the potential. Even excluding 
the complex-mass solutions the perturbative Hamiltonian would appear to be unbounded from below as a 
result of the tachyonic instability. 

^"Here we mean "stable" in the sense of the Ostrogradski instability. The theory contains a tachyon and 
hence it is unstable in this sense. 
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A similar analysis can be performed in the case of open-closed p-adic string theory, a 
simple generalization of ( |1 . 1[ ) which incorporates both the open- and closed-string tachyons. 
This analysis is reported on in appendix C. 

4.3 SFT Near the False Vacuum 



Equation ( |1.2| ) has two constant solutions: (j) = —1 and (/> = 0. The former corresponds 
to the unstable (tachyonic) maximum while the latter is the true vacuum of the theory. 
Writing <j) = —l + 5(j) and linearizing in 6(j) we obtain 

6-^=° [1 + □] = (4.8) 

so that the kinetic function 

F{z) = {I + z) (4.9) 

has only one zero at z = —1: exactly the tachyonic mode which is expected on physical 
grounds. It follows that the initial value problem for ( [4.8D is well-posed with only two 
initial data and the solutions are the tachyonic modes ( p.l4| ) with /i^ = 1. 

4.4 SFT Near the True Vacuum 

We now write cj) = ^ + dcj) and linearize ( |1.2| ) with the result 

[(1 + 0)6"^° - 2] (^0 = (4.10) 

The kinetic function is 

F[z) = e-^'{l + z)-2 (4.11) 
The transcendental equation F{zn) = has solutions 

Zn = -l--Wn{-2ce-') (4.12) 

where Wn are the branches of the Lambert W-function. It is straightforward to check 
that for c = IniS^/i"^) none of the real so that ( [4.10| ) admits infinitely many initial 



data. The solutions are the complex-mass states ( p. IS] ) which, perhaps, can be physically 
interpreted as closed string excitations (at large n equation (4.12) describes a spectrum 



with equally spaced masses-squared). However, such an interpretation should be employed 



with great care since the Hamiltonian for this system is unbounded from below |46]. See 



also [21] for further discussion. 

Notice that, as in subsection we could project out the unwanted complex-mass 
states by taking the contour C in ( |2.3D to enclose only the real and imaginary axes in the 
complex plane. 

5. Nonlinear Equations of Infinite Order 
5.1 The Nonlinear Problem 

A systematic treatment of nonlinear differential equations of infinite order is beyond the 
scope of this paper (indeed, such a treatment is absent even for the case of finite order). 
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Here we discuss some promising approaches and, in the next subsection, we will consider a 
perturbation theory approach to the nonlinear equation of p-adic string theory, rendering 
the initial value problem physically sensible by suitable choice of the contour of integration, 
C. Finally, in subsection 5.3 we will discuss the implications of this choice of C on the full 
nonlinear problem. 

A powerful and efficient method of generating solutions for infinite order nonlinear 
equations (which can be used also in cosmology) is the iterative procedure [30, 31, p^, 39|. 



Iterative techniques are quite common in the mathematical literature on nonlinear integral 
equations and the appeal of such procedures for equations of the form (|l.3| ) is not surprising 
given the relationship detailed in appendix A. For the purposes of characterizing the initial 
value problem, however, such techniques do not seem to be the most appropriate since 
they typically rely on some pre-existing understanding of the asymptotic behaviour of the 
solutions. (One advantage of the perturbative approach described in the next subsection 
is that one can specify the field and its derivatives at t = and allow the system to evolve 
to a unique solution at t = +00.) On the other hand, one might nevertheless argue that 
iterative approaches do still shed light on the initial value problem. This is so because 
for certain stringy models there is numerical evidence that the iterative procedure always 
converges to a unique solution once the sign at infinity and the value of the field at t = 



have been specified |39|. This suggests a class of theories of the form (|l.3| ) with only two 
arbitrary numbers in the solution. 



Another promising approach is the connection to the heat equation y, 27, |40[] which 
may also be used to shed light on the nature of the initial value problem ||2^. However, it is 
not clear if this approach can also be employed for more general kinetic functions which do 
not involve e~"°. In ||2^ solutions in some stringy models were constructed which depend 
only on a finite number of initial data. 

Finally, we should also mention the expansion in eigenfunctions of the kernel which 
was employed in [40| for the case of p-adic string theory and which does not seem to 
have recieved as much attention as the alternative approaches. The expansion in Hermite 
polynomials of |40| seems to provide a very promising approach. Like the perturbative ap- 
proach which we will describe in the next subsection, this expansion is systematic and can, 
in principle, be continued to arbitrarily high accuracy. However, unlike the perturbative 
approach, it does not rely on any a priori assumptions about the initial data. 

At this point we would like to comment upon an apparent discrepency in the literature. 



How does one reconcile the results of [^9[ (also subsections [4.1| and AA) that equations ( |l.l| , 
|1.2| ) admit infinitely many initial data with evidence (coming from the techniques described 
above) that these equations admit only finitely many initial data? These results are not 
necessarily in conflict because many previous studies of equations of the form ( |1.3| ) may 
have (either implicitly or explicitly) placed constraints on the solutions which omit the 
ghost-like states. For example, it was shown in [46| that the Hamiltonian for both p- 
adic string theory and also SFT becomes bounded from below if one demands that the 
tachyon is at the false vacuum as i ^ — cxd. Thus, the ghost-like excitations will not be 
present in any studies which make this assumption. Finally, it was noted in |49] that 
^^The fact that both p-adic string theory and SFT are rendered ghost-free once this auxiliary constraint is 
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perturbation theory projects our the ghost-hke modes in the case of string field theory 
(this is so because in the case of equation ( |1.2D the nonlocal structure of the theory can 
be put into the interaction term by a field re-definition) which explains why path integral 
quantizations do not exhibit any signs of pathology. (See also for details of the 
Feynman rules in string field theory.) 

5.2 Pert ur bat ive Approach 

Consider the full nonlinear equation of p-adic string theory (|1.1| ) with initial conditions 
close to (/> = 1. For simplicity we consider only the homogeneous case = (j){t) (inhomoge- 
nous solutions in p-adic string theory have been considered in | ]68| , |69[). In this case it is 
appropriate to expand the p-adic scalar in perturbation theory as 



oo 



(5.1) 



n=l 



It is straightforward perturb the field equation (1.1) up to second order with the result 



P 



AV2 



P 



6^^U = 

^(2)0 = p{p-l) 



(5.2) 
(5.3) 



In general, for the n-th order perturbation, one will obtain an equation of the form 



P 



P 



5(")0(t) = Jnit) 



(5.4) 



where Jn{t) is constructed from perturbations of order less that n. At each order in 
perturbation theory the generatrix has the form 



/(.)=p^V2 



■P 



(5.5) 



We project out the complex-mass states (discussed in subsection 4J.) by taking the contour 
C to encircle only poles on the real and imaginary axes in the complex plane. Thus we are 
free to write (B.5D in the form 



fis) = j{s)is + V2)is-V2) 



(5.6) 



with 7(s) 



sV2 



p 



at each order in pertur 



/ [s^ — 2] having no zeros in the domain of interest. Therefore, 
:bation theory there are only two initial conditions and hence that 
the full nonperturbative solution also admits only two initial conditions. For the sake of 
illustration we choose 0(0) = 1 + e and (j){0) = where we assume that |e| <C 1. 



imposed means that the ghost-hke modes are not present when one considers the specific physical problem 
of studying brane decay in these theories. However, this resolution seems somewhat unsatisfying since 
one would like to be able to consider such nonlocal theories in a more general setting. The prescription 



for re-defining the pseudo-differential operator discussed in subsection 3.1 provides an auxiliary constraint 



which also projects out the ghost-like states, but is less restrictive than demanding that sit at the unstable 
maximum as t — > — oo. 
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Inspection of (|5.4| ) shows that the full solution will have the form of a sum of expo- 
nentials 

+00 



4>{t) = ^ a„e^ 



similar to what was employed in |13| and |32, |3j]. (In previous applications the terms with 
n < were omitted by the boundary condition (/>(— 00) = 1, here we keep those terms since 
we wish to impose our initial conditions at t = 0.) 

We now proceed to solve the perturbation equations. The linear perturbation { ^.2\j is 
straightforward to solve using ( |2.9D 



= ecosh(^V2t) (5.7) 
which obeys 6^^^(f>{0) = e, 6^^^^{0) = 0. The second order equation ( |5.3| ) takes the form 

p9?/^ -pls^'^U = Mt) (5-8) 



'Mt) 



p{p - 1) [cosh(2\/2t) + 1 



(5.9) 



In Laplace space the source J2is) can be written as 



•^2(5) = jP(p- 1) 



1 



+ 



1 



2 

+ - 



S-2V2 s + s 



(5.10) 



Plugging this into (2.13) and performing the contour integration yields the particular so- 
lution 

S^'^Mt) = 03iy cosh(2^/2i) + '^j^ cosh(^t) - ^ (5.11) 



To (5.11) we are free to add a solution of the homogeneous equation in order to fix the 
initial conditions. The appropriate choice is 



hit) 



p — 1 P — 1 P 
2(p3 - 1) 31np ~ 2 



cosh(\/2t) 



(5.12) 



The full second order solution 6^^'^(p{t) = 6^^Up{t) + 6^^^ obeys 5(2)^(0) = 5(2)^(0) = 
and can be written as 



5(2)<^(t) 



e^(p-l) 
2(p3 _ 1) 



cos h{2V2t) - cosh(\/2t) 



+ 



e^p 



cosh{V2t) - 1 



(5.13) 



In principle this procedure could be continued up to arbitrarily high order in perturbation 
theory. It should be clear the the perturbative method employed in this section could also 
be readily applied to the SFT equation of motion (^]^) or to other nonlinear equations of 
the form ([0| ). 
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5.3 The Implications of Restricting the Contour of Integration 

At the end of subsection we have advocated re-defining the action of the formal pseudo- 



differential operator f{dt)(f>{t) in (2.4) taking C to only enclose the real and imaginary 
axes in the complex plane. We have shown in section ^ that this prescription renders both 
the p-adic string and also SFT ghost-free (at the linearized level) by projecting out the 
unwanted complex-mass states. In subsection |5.2| we have illustrated how this approach 
may be employed order-by-order in perturbation theory so that one could (in principle) 



construct solutions of the full nonlinear problem which depend on only two initial data 



(it should be clear that the same approach will also work for SFT (1^)). However, there is 
nothing inherently perturbative about the definition (|2.4| ) and one could imagine imposing 
this restriction even beyond perturbation theory. Our prescription therefore provides a 
nonperturbative means of re-defining both p-adic string theory and SFT in such a way as 
to evade the Ostrogradski instability.^^ 

However, one may wonder what the implications of this restriction are for the fully 
nonlinear problem. It would be a shame if this re-definition were to project out not only the 
ghosts but also some physically interesting nonperturbative solutions. We now investigate 
this question, showing that our prescription does not affect most of the interesting solutions 



of (1.1, [L^ ) which have been considered in the literature. We will focus primarily on p-adic 
string theory, though it will be clear that similar conclusions will also apply for SFT. 

First, we note that at the nonlinear level, one could think about this restriction on C in 
( |2.4| ) as limiting ourselves only to look for solutions (j){t) within a certain class of functions. 
Specifically, if ^(i) can be written as ( |2.3| ) with C only enclosing the real and imaginary 
axes, then clearly the action of f{dt) is unaltered by our prescription. (This is simply to 
say that if none of the poles {sj} of (f){s) have both Re(si) and Im(si) non-zero then an 



infinite contour C in (2.4) can be deformed to only encircle the real and imagninary axes.) 



This restriction clearly defines a rather large class of functions - it admits any function 
which may be represented as a sum (either discrete or continuous) of exponential modes of 
the form e"^*, e^* with a;, A real. In particular, this restriction does not omit any function 
which has a well-defined Fourier transform (since such functions can be represented as a 
sum of modes e*"^*). Indeed, this restriction is considerably less stringent than demanding 
that the solution (j){t) has a Fourier transform since it also admits a wider class of function 
(such as e*) which do not. 

Let us consider some particular nonperturbative solutions of (|l.l| ) which have been 
constructed in the literature. The anharmonic oscillations of [jl3| (which have been sug- 
gested to represent closed string excitations) can be written as a discrete Fourier series 
and hence this class of solutions is not omitted. Equation (|1.1|) admits a kink-like solution 



p^ , |13| , |37| , 40| which interpolates between the unstable maxima (j) = il at t — > ±00 (for 
odd p). Since this kink solution can be thought of as the zero- frequency limit of the an- 
harmonic oscillations ||l^ it follows that this solution is also untouched by our restriction. 



Similar comments apply for kink-type solutions in SFT and other stringy models |1G]-|19| 



^^The idea of finding some nonperturbative means of restricting SFT to a naive subclass of its solutions 
in order to evade the Ostrogradski instability was suggested in Mm. 
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|3C, 3£]. The rolling solution of (and all solutions which can be constructed 

as summations of the form X]„e"'^*, both in p-adic string theory and in SFT also 
remains unexcluded. Finally, note that equation (^]^) admits the rapidly growing solution 

0(t)=p^expQ^t2) (5.14) 



(see, for example, [40, |69|]). The identity 

" + 00 



implies that this solution can be represented as a sum of exponential functions with real- 
valued argument and hence this solution also survives our restriction. (Rotating t ^ it 
we see that any gaussian lump-type solutions are also unexcluded.) One the other hand, 
clearly our prescription will exclude solutions which involve summations of terms of the 
form ( p. 15 ), such as the infinite-parameter solution of (LI) described by equation (88) of 



fl^ . We see, then, that the restriction which one must impose to render the theories ( |1.1| , 
|1.2| ) ghost-free is not terribly onerous. 

An interesting and important question is whether this prescription preserves the per- 
turbative S-matrix. Though we do not have any rigorous proof, we expect that it does as 
long as one only considers S-matrix elements involving states with real-valued (which, 
in any case, are the physically meaningful states). This agrument is supported by evidence 
from perturbative string field theory. As we have seen in subsection [4.3| , our prescription 
has no effect on SFT if one considers a perturbative expansion about the unstable maxi- 
mum. (This is so because perturbation theory also projects out the unwanted complex-mass 



states [^.) Since perturbative string field theory is believed to be consistent [35| and to 
reproduce the usual string theory amplitudes ||6^ it seems clear that, at least in this par- 
ticular case, our prescription does not unacceptably modify the S-matrix. Of course, this 
argument is far from conclusive, and it is an interesting problem study the implications 
of our prescription on the S-matrix in a more general setting. We leave the complete 
resolution of this issue to future investigations. 



6. Conclusions 

In this paper we have presented a simple and intuitive formalism for studying the initial 
value problem in a wide variety of nonlocal theories. Contrary to naive arguments, dif- 
ferential equations of infinite order do not necessarily admit infinitely many initial data. 
Rather, we have shown that every pole of the propagator contributes two initial data to the 
final solution. Crucially, we have shown that this counting procedure exhausts all possible 
solutions at the linear level. This result has a transparent physical interpretation since 
each pole of the propagator should correspond to a physical excitation and the two initial 
data per physical state are promoted to annihilation/creation operators in the quantum 
theory. 

We have considered, in particular, the dynamical equations of p-adic string theory and 
string field theory arguing that in both cases the initial value problem admits infinitely 
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many initial data. However, both theories may be rendered ghost-free by suitably re- 
defining the action of the pseudo-differential operator (taking the contour in ( |2.4D to enclose 
only the real and imaginary axes) and we have suggested that such a re-definition may be 
analogous to putting a UV cut-off on the theory. However, this procedure lacks firm 
mathematical motivation and may be difficult to interpret physically. Indeed, a skeptic 
might argue that it is tantamount to simply excluding the unwanted solutions by hand. 

On the other hand, we have argued that one might take the contour C as part of the 
definition of the nonlocal theory so that different choices of C yield different theories with 
different mass spectra. In [^] it was suggested that one might evade the Ostrogradski 
instability in string field theory by somehow re-defining the theory to constrain it to a 
naive subclass of its solutions. It was noted there that defining the theory through a 
perturbative expansion (analogous to the analysis in subsection [4. 3D accomplishes just this. 
Our prescription for redefining the action of the pseudo-differential operator accomplishes 
the same thing, though it is non-perturative (there is nothing inherently perturbative about 
( |2.4| ) and, indeed, an analogous definition is used to study nonlinear equations in ||5^). The 
same prescription works for both p-adic string theory and SFT. Moreover, in both cases 
this prescription seems to preserve the nice features of the theory. We have shown in 
subsection that this restriction does not exclude most of the interesting solutions which 
have been constructed previously in the literature. We believe that this is an idea which 
merits further investigation. 



One limitation of our analysis is that we have, for the most part, restricted our attention 
to linear differential equations. Though this certainly omits some interesting problems, it is 
sufficient for many physical applications since a great deal of information can be extracted 
by perturbing about some vacuum solution. We have illustrated this by studying the 
dynamical equation of p-adic string theory up to second order in perturbation theory in 
subsection 5.2. Since, at least in principle, one could construct nonperturbative solutions by 
resumming the perturbation series we expect that a similar counting of initial data should 
apply also for nonlinear equations. However, we leave a detailed study of the solutions of 
nonlinear equations of infinite order to future investigations. 



It is sometimes argued that obtaining any continuum field theory of quantum gravity 
will require abandoning locality |49]. Since presumably any viable nonlocal theory should 
admit only finitely many initial data (at most two if we wish to evade the Ostrogradski 
instability) it is therefore an interesting question as to whether one can construct interacting 
fundamentally nonlocal theories which admit only finitely many initial data (preferably 
without having to resort to re-defining the action of the pseudo-differential operator) . This 
is certainly possible, the model of [EO] provides an explicit example (see also discussions 



in ||47|, ^). We have noted that equations involving fractional operators also seem 
to provide a promising class of theories. We believe that the formalism which we have 
presented in this paper will be useful in finding also more general classes of theories which 
satisfy these criterion. 
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Note Added 



After completing this paper we became aware of work by Lee and Wick |70| and also 
by Cutkosky, Landshoff, Olive and Polkinghorne |7l[] which may have some bearing on 
the our prescription for re-defining the formal pseudo-differential operator. In [^] Lee 
and Wick proposed a modification of electrodynamics which involves finitely many higher 
derivatives. The Lee- Wick model has recently been revived as a possible solution to the 



hierarchy problem |72|. In this construction the propagator for each standard model (SM) 
field contains two poles: one corresponding to the usual SM state and the other corre- 
sponding to an associated Lee- Wick particle. In a modification of the usual contour 
prescription for Feynman diagrams was proposed which preserves the unitarity of the the- 
ory and removes the exponential growth of disturbances associated with the Ostrograski 



instability coming from the Lee- Wick parteners. The modification of |71] seems similar in 
spirit to our suggestion for deforming the contour of integration in the definition of the 
action of the formal pseudo-differential on Laplace space. The prescription employed in 



1 71] can be shown to be equivalent to imposing a future boundary condition that forbids 
outgoing exponentially growing modes and leads to violations of causality. These violations 
of causality are restricted to microscopic scales [73| and the theory is believed to respect 



macroscopic causality and to be free of paradoxes. Though it remains to be seen whether 
similar statements can be made about the prescription which we have advocated, the con- 
sistency of the Lee- Wick model seems to provide a good reason to be optimistic on this 
front. It would be interesting to investigate this issue in detail. 
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APPENDIX A: Convolution Form of the p-adic Equation 



The relationship between equation (|l.l| ) and a certain nonlinear convolution equation has 
previous been noted in the literature |jl3|. Using ( p.4| ) the quantity p^t/^</>(t) on the right- 
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hand-side of can be written as 



1 

1 

27d 



c 



ds (j)(s)e 



St 



C 



1 



+ 00 



1 



y/2Trlnp 
1 



^^/g-(t-t')V(21np) 



g-(t-t')V(21np)g-.(t-t') 



■oo 

+ 00 



ds (i>(s)e 



St' 



c 



^27r Inp 

Then equation (^]^) is equivalent to the following nonlinear Fredholm equation 

i / ~{t-t'f/(2\np)<(.> 



\/27r Inp y_ 



(A-1) 



which gives the convolution form of (See, for example, |13, 37, 4C] for details.) 

The above derivation relied heavily on the gaussian form of the kinetic function in 
equation ( |1 . 1| ) . The action of kinetic operators involving more general functions may also 
be written as a convolution with some kernel. To see this, we write 



+00 

dt' K{t')<p{t + t') = f{dM{t) 

-oo 

If we assume that the integrand of ( |A-2D may be expanded as 

ri=0 



(A-2) 



(A-3) 



then, employing the series expansion for the generatrix (2.2), we see that both sides of 
( A- 2D match for 



/{n)(0) 



+00 



dteK{t) 



(A-4) 



APPENDIX B: Review of Laplace Transforms 

We define a function f{t) and its Laplace transform f{s) by the transformations 

fa+ioo 



1 ra+too 

m = C-'{Hs)} = 17- / e''Ks)ds (B-1) 



oo 



4>{s) = £{(/)(*)} = / e-''(t){t)dt (B-2) 



The ds integration in (B-1) is performed along a vertical line in the complex s-plane and 
a should be chosen sufficiently large that all poles of the integrand are to the left of the 
contour. For f > we can close the contour to the left using an infinite semi-circle so that 

Mt) = — i e'^4>(s)ds for t>0 
2m Jc 
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which gives equation ( |2.3D . On the other hand, for t < the contour should be closed 
to the right and the integration gives zero because the integrand is everywhere analytic 
within the contour of integration. Strictly speaking ( |B-2| ) applies only for Re(s) > a. For 
Re(s) < a the Laplace-space function (/>(s) should be defined by analytic continuation. 
Differentiation acts in Laplace space as 



1 



St 



i=l 



(B-3) 



where 



d, = 5p-^V(0) (B-4) 

and X^j=i is identically zero. Following |58] one might define the action of f{dt)(j){t), with 
f{s) given by (O), as 



/(W(t) = ^ / dse 



St 



n=0 1=1 



(B-5) 



As discussed in section g, in solving (|2.1|) the second term in (|B-5D may be adsorbed into 
the arbitrary coefficients in the solution (/>(s). This is material is discussed, for example, in 



mm 



APPENDIX C: Open-Closed p-adic Strings 



A simple generalization of open p-adic string theory (LI) which incorporates also the closed 



string tachyon is the coupled infinite order system |75, 68 1 

P 



2p 



(C-1) 
(C-2) 



where cp represents the open-string tachyon, ip represents the closed string tachyon and 
<C 1 is related to the string coupling. Let us put the open-string tachyon in its vacuum 
state 4> = and consider the resulting closed-string tachyon dynamics. For = equation 
(|C-2|) becomes 



P 



A^(P-I) 
2p 



^p(p-l)/2-l 



(C-3) 



This equation admits two constant solutions tp = ipf and ip = tpt representing the false 
and true vacuum states respectively. The false vacuum is 



-0/ = 1 + 



A2 



2p{p + 1 

Writing ip = ipf + 6ip and linearizing ( |C-3| ) in 6ip gives 



p' + y 



4p 



61^ = 



(C-4) 



(C-5) 
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to leading order in A^. The kinetic function 

F{z) = v~"'^ - ( + A 
has infinitely many zeros at z = z„ where 



2 , x2 P'^+P-2 
4p 



Inp 



In 



V 



Ap 



+ 



Snin 
Inp 



(C-6) 



(C-7) 



and n = 0, ±1, • • • 

We consider now fluctuations about the true vacuum ipf For p > 2 (the case p = 2 is 
treated below) we have 

V'. = (C-8) 

and writing ip = + dip yields 

p-°/45^ = (C-9) 

The kinetic function F{z) = p~^l^ has no zeros and hence there is no solution. For p = 2 
the true vacuum is shifted by the interaction 



Ipv 



(C-IO) 



however, writing -0 = — + 5^ yields again F{z) = p~^^^. 

One might instead imagine writing (j) = \ + 5(f) and %[) = 1 + 5-0 in which case the system 
of equations ( |C-2 , C-1 ) gives 



-□/2 
p ' — p 



-□/4 2 
p ' — p 



5^ 



'a <^ 

Ap 



(C-11) 
(C-12) 



For A^ = the study of this system is very similar to the analysis in section ^. For finite 
A^ solutions are also possible and have been described in [ p8| . 
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